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We consider the problem of searching a d-dimensional lattice of A'^ sites for a single marked 
location. We present a Hamiltonian that solves this problem in time of order \fN for d > 2 and 
of order \/N\ogN in the critical dimension d = 2. This improves upon the performance of our 
previous quantum walk search algorithm (which has a critical dimension of d = 4) , and matches the 
performance of a corresponding discrete-time quantum walk algorithm. The improvement uses a 
lattice version of the Dirac Hamiltonian, and thus requires the introduction of spin (or coin) degrees 
of freedom. 



I. INTRODUCTION 

Quantum mechanical computers can solve certain 
problems asymptotically faster than classical computers. 
One of the major advantages of quantum computation 
comes from a fast algorithm for the problem of finding a 
marked item among N items. Whereas a classical com- 
puter requires Q{N) steps to solve this problem, Grover 
showed that a quantum computer can solve it using only 
0{\/N) steps Xl, which is optimal |2|- 

To apply Grover's algorithm, it must be possible to 
quickly perform a reflection about a superposition of all 
possible items. However, this may not be feasible if the 
items are distributed in space and the algorithm is re- 
stricted to access them by local moves. For example, if 
the items are arranged on a one-dimensional line, simply 
traveling from one end of the line to the other requires 
N moves, and a straightforward argument shows that no 
local algorithm, classical or quantum, can find a marked 
item in less time than Vl{N). But for other geometries, 
such as higher dimensional lattices, a quantum algorithm 
can conceivably achieve a speedup over the classical com- 
plexity of e(iV). 

Recently, there has been considerable progress in un- 
derstanding the spatial search problem for quantum com- 
puters. Aaronson and Ambainis gave an algorithm that 
finds a marked item in the optimal time 0{\fN) for a lat- 



tice in d > 2 dimensions, and in time 0( V N log TV) for 
a two-dimensional lattice Their algorithm is based 
on a carefully optimized recursive search of subcubes, 
which raises the question of whether a simpler algorithm 
could solve the problem just as quickly (or perhaps even 
faster in two dimensions). In particular, it is interest- 
ing to consider quantum walk algorithms, which only use 
local, time-independent dynamics. Two distinct kinds 
of quantum walk algorithms have been considered. In 
the continuous-time quantum walk the algorithm is 
described by a time-independent Hamiltonian connect- 
ing adjacent sites. In the discrete-time quantum walk 
UTIE Q J the algorithm consists of repeated application 



of a fixed local unitary transformation. 

In , we considered a continuous-time quantum walk 
algorithm for the spatial search problem using no addi- 
tional memory beyond the present location of the walker. 
We showed that this algorithm can find a single marked 
site in time 0{^/N) for dimensions d > 4 and in time 
0{\/N log N) in four dimensions.^ We also showed that 
this algorithm fails to provide an interesting speedup for 
dimensions d < 4. Subsequently, Ambainis, Kempe, and 
Rivosh found a discrete-time quantum walk algorithm 
that works in lower dimensions ■ This algorithm runs 
in time 0{Vn) for d > 2 and in time 0{Vn log N) in 
two dimensions. Because a discrete-time quantum walk 
cannot be defined on a state space consisting only of the 
vertices of a graph the algorithm of [l^ necessar- 
ily uses additional memory (sometimes referred to as a 
"coin" in analogy to classical random walks). In this 
paper, we consider a continuous-time quantum walk us- 
ing additional memory, and we show that it achieves the 
same running times as the discrete-time algorithm. 

Through the analysis of i^], the failure of the 
continuous- time quantum walk algorithm for d < 4 can 
be viewed as a consequence of a quadratic dispersion re- 
lation: states with small momentum have energy pro- 
portional to their momentum squared. If the dispersion 
were linear instead of quadratic, so that states with small 
momentum had energy proportional to their momentum, 
one might expect the algorithm to work whenever d > 2. 
A natural way to achieve linear dispersion is to employ 
the massless Dirac equation. Indeed, using an appropri- 
ate lattice version of the Dirac Hamiltonian, we find a fast 
algorithm for spatial search in d > 2. Because the Dirac 
particle necessarily possesses spin degrees of freedom, the 
resulting algorithm must have additional memory beyond 
the present location. Thus, although a continuous-time 
quantum walk can be defined without additional mem- 
ory, we find that the additional degrees of freedom can 
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^ The run time 0{\/~N log N) in d = 4 can be achieved using am- 
plitude amplification Using only classical repetition of the 
quantum walk, the algorithm requires 0(v']V log^^^ A'^) steps. 
The same remark applies to the discrete-time quantum walk al- 
gorithm in d = 2. 
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improve the algorithm's performance. 

II. HAMILTONIANS FOR SPATIAL SEARCH 

In the Hamiltonian formulation of the spatial search 
problem, our goal is to write down a local Hamiltonian 
that will quickly transform a simple initial state, such as 
the uniform superposition over all lattice sites 



(1) 



to a state with substantial overlap on the marked state 
\w). In M, we considered the Hamiltonian 



To write down the Dirac equation in d dimensions, 
we need d + 1 anticommuting operators. The minimal 
representation of the algebra (0) uses 2 T'^/^l -dimensional 
matrices, and hence there are 2'^'^/^^ spin components. 

Now consider a lattice version of the massless Dirac 
Hamiltonian, equation ((Sj with m = 0. The continuum 
operator pj can be discretely approximated as 



(7) 



where ej is a unit vector in the j direction. However, 
as we will see later, it turns out that simply taking the 
free Hamiltonian @ using the lattice approximation (0 
is insufficient. Instead, we will take^ 



= -jL- \w){w\ 



(2) 



Here the second term identifies the marked location, 7 
is an adjustable parameter, and L is the Laplacian of an 
A'^-site square lattice in d dimensions, periodic in each 
direction with period N^^'^. L has matrix elements 



{x'\L\x) 



1 X adjacent to x' 
—2d x ~ x' 
otherwise. 



(3) 



It is called the Laplacian because it is a discrete approx- 
imation to the continuum operator V^. 

Since the free Hamiltonian — 7L is translationally in- 
variant, its eigenstates are the momentum eigenstates 



\k) 



1 



E 



(4) 



where kj = 

N^/'^ odd, and m 



withmj = 0,±1,. ..,±i(iVi/''-I) for 
.J = 0,±l,...,±i(iVi/rf - 2),+i7Vi/'^ 
for iV^/'' even. The energy of the state \k) is £L{k) — 
27(1 — X]j=iCOsfcj). Thus, for small |fc|, £L{k) « 'fk'^. 
This quadratic dispersion relation ultimately gives rise to 
the critical dimension d = 4 for the algorithm of . 

To achieve linear dispersion, we can replace — by 
the massless Dirac Hamiltonian. In the general case of 
mass m, this Hamiltonian has the form 



Dirac 



(5) 



(8) 



where both uj and 7 are adjustable parameters. For a 
Hamiltonian with spin degrees of freedom, translation 
invariance shows that the eigenstates have the form jry, fc), 
where |ry) is a (momentum-dependent) spin state. For 
©, we find states with energies 



£{k) = ±y/uj^s^{k)+j^c^{k).^ 



(9) 



where 



s^{k) ^"^sm"^ kj , c(fc) 2^(1 - cos fcj) . (10) 



For small momenta, we have £{k) w ±w|fc|, which leads 
to a better search algorithm in low dimensions. 

The full algorithm is as follows. We begin in the state 
|ry, s), where \ri) is any spin state and \s) is the uniform 
superposition We then evolve with the Hamiltonian 



H = Ho-P\w){w\ 



(11) 



with parameters ^,7 to be determined in the analysis 
below, for a time T also determined below. The goal is 
to choose the parameters uj and 7 so that for some T as 
small as possible, the spatial component of the evolved 
state has a substantial overlap on \w). 



where the operators aj and /3 act on spin degrees of free- 



dom, and p 



' da; 



is the momentum operator. If the 



spin operators and /3 satisfy the anticommutation re- 
lations 



{a„/3} = 0, /3^-l, (6) 
m , 

0, i?Diiac has linear 



then one finds -ff^irac ~ 1^1^+ ^ required for a rel- 
ativistic particle. Thus for m 
dispersion, i?Diiac = ±bl- 



III. ANALYSIS OF THE ALGORITHM 

To analyze the algorithm, we would like to determine 
the spectrum of H using our knowledge of the spectrum 
of Hq. We do this using the same techniques we applied 



^ This choice is closely related to a standar d re medy for the fermion 
doubling problem in lattice field theory Il3l p. 27]. 
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to the Hamiltonian in An eigenvector of H, de- 
noted l^pa), with eigenvalue Ea, satisfies 

H\^Pa) = {Ho - f3\w){w\)\^ba) = Eal^Pa) ■ (12) 

Defining 

Hi'a) = /R^|0a) (13) 

where {(fia) is a normalized spin state, and \/Ra > by 
choice of phases, we can rewrite (|12|) as 

(Ho-EaMa) =^/R^ma,w). (14) 

Assuming Hq — Ea is nonsingular, we can write the eigen- 
state of H as 

m = J^^^ma,w). (15) 

Consistency with (|13() then gives the eigenvalue condition 




FIG. 1: Critical values of (0^,7) for various dimensions. From 
rightmost curve to leftmost curve, d = 2, 3, 4, 5. 



where 



|0a) = F{Ea)P\cPa) 

F{E) - {w\ — 

-no 



1 



(16) 



(17) 



operates on the spin degree of freedom. To find eigenval- 
ues oi H , we must look for values of E such that the spin 
operator F{E)P has an eigenvector of eigenvalue 1. 

In addition to finding eigenvalues of H, we need some 
facts about its eigenvectors. The normalization condition 
on \-)pa) gives 

{^^,w\p—^-—^/3\(l)a,w) (18) 



R„ 



{Ho ~ Ea)^ 
{^a\PF'{Ea)ma). 



(19) 



We also need the overlap of {ipa) with eigenvectors of Hq- 
From (|15|l we have 



{£\^a) 



'Ra 



£-Ea 



{£ma,w) 



where \£) is an eigenvector of Ho with eigenvalue £. 
For the free Hamiltonian (|SJ|, we find 



F{E)p 



Ho+E 



'Hi 



7c(fc)+/3E 

£;2 



N £{kY 



^~^ + U{E)+pEV{E), 



(20) 

(21) 
(22) 

(23) 



where in H22|l we have canceled terms that are odd in fc, 
and 



7 c{k) 



N ^ £(fc)2 



E^ 



N V £{kV 



E^ 



(24) 
(25) 



If \E\ ^ £{k) for all k 0, then we can Taylor expand 
U{E) and V{E) in powers of E. We will need only the 
leading order terms U{0) and V{0). For large N, we have 



U{0) 



-fc{k)d''k 



{2^YJ_^ £{ky 



(26) 



which is a convergent integral regardless of d. For d > 2 
and N large, we can also write V{0) as a convergent 
integral, 



V{0) 



1 



{2nr £{kf 



(27) 



In d = 2, this integral is logarithmically infrared diver- 
gent, and instead we find 



no)- 



1 



47rw 



■ln7V + 0(l). 



(28) 



Now suppose we choose lj and 7 such that C/(0) = 1. 
In this case, we are simply looking for a zero eigenvalue 
of l3{—j^ + EV{Q j). We find such eigenvalues with 



E±^± 



(29) 



which indeed satisfy the condition \E±\ ^ £{k) for all 
/c 7^ 0. These eigenvalues are degenerate in the spin 
space, i.e., any state \4>±) provides an eigenvector with 
the same eigenvalue. 

The condition U{0) = 1 can be satisfied by choosing 
u{lu/j) = 7, where u{lu/j) = "/U{0) is a function only 
of uj/'^. Figure n shows the critical curve in the {uj,'y) 
plane for d — 2 through 5. For any uj with < u; < w*, 
where uj* is some dimension-dependent threshold value, 
there are two values of 7 such that U{Q) = 1. Note that 
with = 0, we recover the results of M.. Also, with 
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7 = 0, no solution of U{0) = 1 exists, so it was essential 
to include this additional term (and the ability to fine 
tune its coefficient). 

Having found the relevant eigenvalues, we need to de- 
termine the corresponding eigenstates. Using H19|l we 
find 



i?4 



f 



NEl 



y(o) « 2V{0) , 



and using 12U|) we find 

(7?,S|V'±) = - 



E±VN 



1 



where we have chosen the eigenstate of H with \(i>±) 
/3|?7). Therefore we have 

1 



In high dimensions, our algorithm is very wasteful in 
terms of the number of spin degrees of freedom: it uses 
a 2 T'^/^l -dimensional spin space, whereas shows that 
no spin degrees of freedom are required at all for d > 4. 
In comparison, the discrete-time quantum walk search 
algorithm in [l^l uses 2d extra degrees of freedom. The 
Dirac particle in d dimensions cannot be represented with 

(30) fewer than 2'^'^/'^'^ de grees of freedom, but a continuous- 
time search algorithm with only d+1 degrees of freedom 
can arise from reproducing the Dirac algebra © only on 
a subspace. If the operators aj and j3 satisfy 

(31) afcll'?) = 25,- fcl'?) , {a,.P}\ri) = , /3h) = (34) 

for some spin state then the algorithm will work 
starting from the state jyy, s). The condition H34|) is suffi- 
cient to give |?7, k) — £{ky'\rj, k). The previous analy- 
sis then shows that 



|r;,s)«-^(|V'-)-|V^+)), 

and choosing T ~ tt/{2\E±\) produces the state 

1 



-iHT 



(32) 



(33) 



which has an overlap on \ri,w) of ^J2R±. 

For d > 2, we have shown that there is a T = 0{^/N) 
that gives a probability 0(1) of finding w. For d = 
2, there is a T = 0{y'N log N) that gives an ampli- 
tude 0{l/\/\ogN), so that amplitude amplification 
can be used to find w with a probability 0(1) in time 
0{VN\ogN). 



IV. DISCUSSION 

We have described a continuous-time quantum walk 
algorithm for the spatial search problem. Using Dirac's 
insight of introducing spin to take the square root in a 
relativistic dispersion relation, we have found a Hamil- 
tonian that locates a single marked item in the optimal 
time of 0{VN) above the critical dimension (d > 2), and 
that runs in time 0{^/N log N) in d = 2. 

This algorithm is closely related to the discrete-time 
quantum walk search algorithm of [To| . Very similar tech- 
niques to the ones we have used in this paper can also be 
applied to discrete-time quantum walks 'l4'|. This analy- 
sis for the algorithm of 10] closely parallels the analysis 
above, which highlights the similarity between the two 
kinds of algorithms. However, there are a few important 
differences. The continuous-time algorithm requires fine 
tuning the parameters uj and 7, whereas there is (appar- 
ently) no equivalent fine tuning in the discrete-time algo- 
rithm. Also, the discrete-time algorithm has noticeably 
different behavior depending on whether A^^/'' is odd or 
even, a difference that is not seen in the continuous-time 
algorithm. In short, although the essential infrared fea- 
tures of the two kinds of algorithms are identical, their 
detailed behaviors differ. 



'Ra 



Hq — Ea 



(35) 



is an eigenstate of H with eigenvalue Ea provided 
-ivk + U{Ea) + EaV{Ea) = 1, whcrc U{E) and V{E) 
are as defined in equations H24|) and H25|l . The rest 
of the analysis with two states jf/'i) follows exactly as 
before. Finally, we see that (|34|l can be satisfied in a 
(d+ l)-dimensional spin space with basis |0), |1), . . . , |d), 
since in that case we can choose aj ~ |0)(j| + b)(0|. 



/3 = 2|0)(0|-/, and h) = 10). 

Unlike the algorithm of [g , the algorithm of this paper 
cannot be turned into an adiabatic algorithm. With the 
Hamiltonian by starting in the state \s) and lower- 
ing the parameter 7 from a large value to zero sufficiently 
slowly, the adiabatic theorem guarantees that the system 
will remain near its ground state, ending up close to the 
state \w). In d > 4, this can be done in time 0{^fN)^ 
and in d = 4, it can be done in time 0{-/N\o^^'^ N). 
However, in the algorithm of the present paper, states 
with k — Q are not the ground state of the free Hamil- 
tonian ©; these states have zero energy, but this is in 
the middle of the spectrum. Although the adiabatic the- 
orem applies to any eigenstate, not just the ground state, 
states near the middle of the spectrum of Hll|) with w,7 
small have very little overlap on jw), so that even per- 
fectly adiabatic evolution produces a state far from the 
desired one. 

Finally, we note that the actual complexity of the spa- 
tial search problem in d = 2 is still an open question. A 
gap of log A'^ remains between the best known algorithm 
and the lower bound of 0]. It would be interesting to 
improve the algorithm further or to show that no such 
improvement is possible. 
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